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Abstract 



We compute lattice renormalisation constants of one-link quark operators (i.e. 
operators with one covariant derivative) for overlap fermions and Liischer-Weisz 
gauge action in one-loop perturbation theory. Among others, such operators enter 
the calculation of moments of polarised and unpolarised hadron structure functions. 
Results are given for (3 = 8.45, (3 = 8.0 and mass parameter p = 1.4, which are com- 
monly used in numerical simulations. We apply mean field (tadpole) improvement 
to our results. 
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1 Introduction 



In a recent publication pQ we have computed lattice renormalisation constants of local 
bilinear quark operators for overlap fermions and improved gauge actions in one-loop 
perturbation theory. Among the actions we considered were the Symanzik, Liischer- 
Weisz, Iwasaki and DBW2 gauge actions. The results were given for a variety of p 
parameters. Furthermore, we showed how to apply mean field (tadpole) improvement 
to overlap fermions. In this letter we shall extend our work to one-link bilinear quark 
operators. Operators of this kind enter, for example, the calculation of moments of 
polarised and unpolarised hadron structure functions. The present calculations are much 
more involved than the previous ones, so that we shall restrict ourselves to the Luscher- 
Weisz action, and to parameters actually being used in numerical calculations. 

The integral part of the overlap fermion action j2J El H] 



where Dy/ is the Wilson-Dirac operator, and p is a real parameter corresponding to a 
negative mass term. At tree level < p < 2r, where r is the Wilson parameter. We take 
r = 1 and consider massless quarks. 

Numerical simulations of overlap fermions are significantly more costly than simula- 
tions of Wilson fermions. The cost of overlap fermions is largely determined by the condi- 
tion number of X'X. This number is greatly reduced for improved gauge field actions [5]. 
For example, for the tadpole improved Luscher-Weisz action we found a reduction factor 
of > 3 compared to the Wilson gauge field action The reason is that the Luscher-Weisz 
action suppresses unphysical zero modes, sometimes called dislocations |7]. A reduction 
of the number of small modes of X*X appears also to result in an improvement of the 
locality of the overlap operator |3] . 

We consider the tadpole improved Luscher-Weisz action (HI El HH| 




(1) 



m being the mass of the quark, is the Neuberger-Dirac operator 
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plaquctte 



rectangle 



(3) 




2 



where £7pi aque tte is the standard plaquette, {/rectangle denotes the loop of link matrices 
around the 1x2 rectangle, and f/ pa raieiiogram denotes the loop along the edges of the three- 
dimensional cube. It is required that cq + 8ci + 8C3 = 1 in the limit g — > 0, in order to 
ensure the correct continuum limit. We define 



= -W (4) 
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The remaining parameters are jlOj : 



ci (1 + 0.4805 a) c 3 0.03325 a 1 ( c x c 3 \ 

— = 1 + 8 — + — ' ( 5 ) 



c 20 ul c ul c V c o ' c o 

where 



M ° = ' 3 Tr ) ' a = ~3 06839 ' ^ 



The final results cannot be expressed in analytic form (as a function of (3 and p) 
anymore. We therefore have to make a choice. Here we consider two couplings, (3 = 8.45 
and 8.0, at which we run Monte Carlo simulations at present jHlE]- The corresponding 
values of c\ and c 3 are [12] : 



p 


Cl 


c 3 


r /a 


8.45 


-0.154? 


546 


-0.0134070 


5.29(7) 


8.0 


-0.169? 


305 


-0.0163414 


3.69(4) 



In (jZJ) we also quote the corresponding force parameters r /a, as given in Assuming 
that r = 0.5 fm, they translate into a lattice spacing of a = 0.095 fm at (3 — 8.45 and 
a = 0.136 fm at /3 = 8.0. The mass parameter was chosen to be p = 1.4. This appeared 
to be a fair compromise between optimising the condition number of X^X as well as the 
locality properties of D N [T5] . 

The paper is organised as follows. In Section 2 we give a brief outline of our calculations 
and present results for the renormalisation constants in one-loop perturbation theory. In 
Section 3 we tadpole improve our results, and in Section 4 we give our conclusions. 



2 Outline of the calculation and one-loop results 

The Feynman rules specific for overlap fermions jTJJ are collected in [T|, while the 
gluon-operator and the gluon-gluon-operator vertices (needed for the cockscomb and op- 
erator tadpole diagrams) are independent of the fermion action and can be found in [16 . 
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We consider general covariant gauges, specified by the gauge parameter £. The Landau 
gauge corresponds to £ = 1, while the Feynman gauge corresponds to £ = 0. In lattice 
momentum space the gluon propagator D pu (k) is given by the set of linear equations 



E 

p 

where 



Gfj,p(k) — - — - k p k p 



D pv {k) = 5^ , (8) 



G^(k) — kpk v + 'y ^ (kpd^v k p k p 5 pi ^j dp P (9) 
p 

and 

d pv = (1 - <W [C - Ci a 2 P - C 2 a 2 (^ + fc 2 )j , fc M = - sin ^ , P = ]T fc 2 . ( 10 ) 

The coefficients {Cj} are related to the coefficients {q} of the improved action by 

(70 = 00 + 801 + 803, Ci = c 3 , C 2 = c 1 -c 3 . (II) 

The calculations are done analytically as far as this is possible using Mathematica. Part 
of the numerical results have been checked by an independent routine. 

The bare lattice operators 0(a) are, in general, divergent as a — > 0. We define finite 
renormalised operators by 

O s (fi) = Z$(a^)0(a), (12) 

where S denotes the renormalisation scheme. We have assumed that the operators do 
not mix under renormalisation, which is the case for the operators considered in this 
letter. The renormalisation constants Zq are often determined in the MOM scheme first 
from the gauge fixed quark propagator Sn and the amputated Green function A of the 
operator O: 

Z^ OM (a^)S N \ p2=ti2 = S tree , (13) 

^MOMf a \ 

° MnM , A I o 2 = A« cc + other Dirac structures . (14) 
z i? (a, A*) p M 



(Note that Z^ = l/Z 2 .) The renormalisation constants can be converted to the MS 
scheme, 

7~MS(„ ,A — yMS,MOM yMOM f ,,\ 

z^, (a,//j— z, z^ (a, /ij, 

15 

Z^ 5 (a, /i) = Z% s ' MOM Zg OM (a, /i) , 
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where z¥ ' MOM , z^ s ' MOM are calculable in continuum perturbation theory, and there- 
fore are independent of the particular choice of lattice gauge and fermion actions. 



In the wave function renormalisation constants where found to be 

ZM° M (a, /i) = 1 - ^ [2(1 - logM + 4.79201 £ + fc] (16) 
in the MOM scheme, and 

Z$^(o, /i) = 1 - f-4 [2(1 - log(a/x) + 3.79201 £ + &e + 1] (17) 

Y l07T Z 



in the MS* scheme, with Cp = 4/3 and 



p 




8.45 


-17.429 


8.0 


-17.054 



We consider the following one-link operators 

= Dvi/>(x) , (19) 

Ojl = |^)7„75 B„ V(^) , (20) 

<— > — > < — 

where D V =D V — D v is the (symmetric) lattice covariant derivative. While in our previous 
work lj, which involved local bilinear quark operators, we only had to deal with the vertex 
diagram shown on the left-hand side of Fig.^ we now obtain contributions from additional 
diagrams: the operator tadpole and the cockscomb diagrams shown on the right-hand side 
of Fig. □ 




Figure V. The one-loop lattice Feynman diagrams contributing to the amputated Green 
function. From left to right: vertex, operator tadpole and cockscomb diagrams. 
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Action 




b 2 


bs 


&4 


P = 8.45 


- 5.6115 


-3.8336 


2.7793 


0.3446 


(3 = 8.0 


- 5.2883 


-3.7636 


2.7310 


0.3331 


Plaquette 


-10.6882 


-4.7977 


3.4612 


0.5267 



Table 1: The coefficients for the tadpole improved Liischer-Weisz action at (3 — 8.45 
and 8.0, as well as for plaquette action. 



The amputated Green function of the operator 0^ u [eq. (fTTIJ) ] turns out to be 



A^(a,p) = f y ( j)u + 



9 2 C, 
16tt 2 



l + n\og(a 2 p 2 ) 



4.29201 C + 



+ 



^log(oV) + ^ + 62 



+ b 4 buvlvVv + 



1+t 



lvV» + 

PfiPu 
p2 



--log(aV) 



<W (21) 



where p is the external quark momentum, and the coefficients {bi} are given in Tabled 
for the tadpole improved Liischer-Weisz action and, for comparison, for the plaquette 
action (with c\ = C3 = 0) as well. The latter numbers are independent of (3. The Green 
function h b {a,p) of the operator O h [eq. (|2~U|l ] is obtained by multiplying the right-hand 
side of (|21|) by 75 from the right. The coefficients turn out to be identical to {h}, 
as is expected for overlap fermions. Thus, 0^ v and O 5 have the same renormalisation 
constants. In the following we may therefore restrict ourselves to the operator O 



/IV • 



It has been checked numerically that the gauge dependent part of (|2Tj) is universal (i.e. 
independent of the lattice gauge and fermion action), in accordance with the arguments 
presented in [T]. 

Under the hypercubic group H{4) the 16 operators of type (|T9~|) fall into the following 
four irreducible representations jTTj: 



.(6) 



a 



= -(0i4 + a 



T 



(3) . 



o = 



T 



r- 



(i) . 
1 

(6) . 



o 



t'2c 



o = 



41, 



044 - - (Oil + 22 + O33) , 

011 + 22 + 033 + 044 , 
014 - 041 • 



(22) 

(23) 

(24) 
(25) 



(We have given one example operator in each representation. A complete basis for each 
representation can be found in [TV| . ) The operators and are widely used in 
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numerical simulations ^HIEHEIE]- They correspond to the first moment of the parton 
distribution. The operators (|24J1 and (|2*3jl represent higher twist contributions in the 
operator product expansion, and so are not used as much as operators in the first two 
representations. For completeness we give results for all four representations, so that the 
renormalisation factors for all operators of the form (fTH|) will be known. We denote the 
corresponding amputated Green functions by A V2a , A V2b , A V2c and A V2d . From (|2*T|) we read 
off 



A 



A 



''2b 



A 



t'2,7 



A 



''I'd 



- ( 7lP4 + 74^1 



1 + 



9 2 C f 



16tt 2 



(£ + -) log(oV)- 3.79201 £ + b V2a 



9 2 C F 

167T 2 
74P4 - 



9 2 C F 



P 



7iPt 



t=l 



1 + 



9 2 C F 
16vr 2 

3 



iog(oV)- 3.79201 e+ 6. 



''26 



16vr 2 



9 2 C F 
+ - — - 

167T 2 
7ip 4 - 74^1 



[(£-1) log(aV)- 4.79201 £ + 6„ 2c ] 



1 + 



9 2 C F 
16tt 2 



[(C - 1) log(a 2 p 2 ) - 4.79201 £ + 6„ 2d ] 



(26) 



(27) 

(28) 
(29) 



with 



61+62 



^26 



6l + 6 2 + 64 



6 1 + 6 2 + 46 3 + 6 4 --, b VM = b x -b 2 . (30) 



It is worth pointing out that with Wilson or clover fermions the Green functions A V2c 
and A V2d both show perturbative mixing of 0(g 2 /a) with local operators. With overlap 
fermions these 0(1 /a) terms are completely absent, showing once again that overlap 
fermions behave much more like continuum fermions when mixing is a possibility. 

Using (fPfj) and (|16|h we obtain the renormalisation constants in the MOM scheme: 



yMOM . 

V2a,V 2 b 

yMOM . 
V 2 c,V2d 



a, jj) 



a, n) 



9 2 Ci 



16tt 2 

F 



16 



9 2 C 



167T 2 



3 



log(a/i) + £ + 
+ 6 S ] . 



+ 6, 



(31) 
(32) 



As already mentioned, the conversion factors Z^ 18 ^ ^ v are universal JTB]. They are 
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given by 



7 MS,MOM 

J V2a,V2b 



7 MS,MOM 



7 MS,MOM 
J v 2 d 



9 2 C F 



IQtt 2 

F 



9 2 C 



16tt 2 



9 ^ 



1 . 



In the MS scheme we then find 



7 MS 

J V2a,V'2b 



a,/i) 



9 2 C F 



167T 2 
F 



16 



40 



9 2 C 



16vr 2 

F 



bg(a/i) + — + b V2a:V2b + 6 S 



S 2 C 



16tt 2 



(33) 

(34) 
(35) 

(36) 
(37) 
(38) 



3 Tadpole improved results 



A detailed discussion of mean field - or tadpole - improvement for overlap fermions and 
extended gauge actions has been given in pQ. Here we will briefly recall the basic idea, 
before presenting our results. 



Tadpole improved renormalisation constants are defined by 

7TI _ 7 MF ( Zo \ 
^O-^O I 7 MF I 

y^O J pert 



(39) 



where Zq F is the mean field approximation of Zq, while the right-hand factor is computed 
in perturbation theory. For overlap fermions (with r = 1), and operators with no covariant 
derivatives, we have 



7MF 



1 — fin 

pu 



(40) 



p-4(l-«o)' 

In our case no = 1- It is required that p > 4(1 — uq), which is fulfilled here (see Table 0). 



To compute the right-hand factor in (|39|h we have to remove the tadpole contribu- 
tions from the perturbative expressions of Zq first. This is achieved if we re-express the 
perturbative series in terms of tadpole improved coefficients: 



c _ 4 ^0 w_ 

2 9 ' 2 

9ti 9 2 9ti 9 
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p 


UTI 




8.45 


0.543338 7T 2 


0.65176 


8.0 


0.515069 7T 2 


0.62107 



Table 2: The coefficient k^ 1 and the average plaquette Uq at (3 = 8.45 and 8.0. 

This does not fix all parameters, but leaves us with some freedom of choice. The simplest 
choice is to define 

a 2 

9ti = —-, c =c , c t =u Ci, 2 = 1,3. (42) 





u 



With this choice 

C " = Co + 8cr + 8cf, C^ugd, C^ = ulC 2 . (43) 

(Note that Cq 1 ^ 1. However, C^f 7 — > 1 in the continuum limit.) This means that we 
have to replace every g 2 by g\j and every c\ and c 3 by of 7 and c 7 / 7 , respectively, while 
keeping c unchanged. The effect of introducing tadpole improved coefficients (JHJ) is that 
the rescaled gluon propagator remains of the same form as we change uq, thus ensuring 
fast convergence. 

To compute Zq F perturbatively, we need to know the perturbative expansion of u 
to one- loop order (201 EI] . We write 

«o = l-%%e- (44) 

167T Z 

In we have computed k FI for the Liischer-Weisz action with coefficients Cq 1 , Cf 1 
and Cvf 7 . The numbers are given in Tabled for our two values of (3, together with the 
'measured' values of Uq. Expanding (|47IJ) then gives 



tMf _ 1 , Qti^f 4 TI 
16tt 2 p u 



« rt = l + ^£f-e- (45) 



Let us now rewrite the one-loop renormalisation constants of Section |2l as 



C F g 

V2a,V2b x in n 



2 



16 

3a; 



log(a/i) +B V2atV2b (p, C) 



(46) 



Cf g 2 

Zv 2c ,V2d = 1 1« 2 Bv2c,V2d(Pl C) ■ (47) 
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Dividing (fUjj) and (}4Tj) by (|43j) and inserting (f4*0"j). we obtain mean field/tadpole improved 
renormalisation constants: 



yTi 



7TI 



P 



p- 4(1 -w ) 



1 - 



1 



16 



167T 2 



J^ti 1o sM + BH jV2b 



p- 4(1 -wo 

where we have introduced the abbreviated notation 

B TI = B{p,C TI 



-k TI 



P 



(48) 
(49) 

(50) 



The coefficients -B(p, C T/ ) are the analogue of B(p, C), with C , Ci and C 2 being replaced 
by Cjf 7 , Cj 1 and Cvf 7 , respectively. In fl4*%j) and (J4T?)) only the gluon propagator has been 
tadpole improved. 



To tadpole improved the fermion propagator as well, we must replace p by [T] 

p- 4(1 -m ) 



WO 



(51) 



in the right-hand perturbative factor of (}3*9"j) . This defines 'fully tadpole improved' renor- 
malisation constants 

16 



7 FTI 

J V2a,V2b 



7 FT I 

J 1>2c,V2d 



p- 4(1 -u ) 
P 



Oti^i 
16tt 2 



3 C " 



log(a/i) + B 



FTI 

V2a,V 2 b 



with 



p- 4(1 -M ) 

B FTI = B(p TI ,C TI ) 



9rfiF tjFTI 
lQ n 2 v 2c ,v 2d 



+ ^TI k u 



TI 



(52) 
(53) 

(54) 



In Table 0] we present our final results and compare tadpole improved and unimproved 
renormalisation constants. We see that the improved coefficients B are rather small in the 
case of the operators t>2 a and V2b, much smaller than for Wilson and clover fermions fU\ . 
which raises hope that the perturbative series converges rapidly. This furthermore means 
that the dominant contribution to the renormalisation constants is given by the mean 
field factor (jUJ). 



4 Summary 

We have computed the renormalisation constants of one-link quark operators for overlap 
fermions and tadpole improved Liischer-Weisz action for two values of the coupling, f3 = 
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Operator 


P 


B 




T>TJ 




n FT J 


/ r J- i i ivi id 




8.45 


-22.430 


1.315 


0.502 


1.393 


-0.077 


1.411 


V2b 


8.45 


-22.085 


1.311 


0.793 


1.384 


0.230 


1.401 


V2c 


8.45 


-18.079 


1.254 


2.985 


1.318 


1.829 


1.353 


V2d 


8.45 


-19.207 


1.270 


2.303 


1.338 


1.369 


1.367 


V2a 


8.0 


-22.036 


1.310 


0.603 


1.390 


-0.108 


1.412 


V2b 


8.0 


-21.703 


1.305 


0.892 


1.381 


0.199 


1.402 


V2c 


8.0 


-17.890 


1.252 


3.038 


1.316 


1.643 


1.358 


V2d 


8.0 


-18.954 


1.266 


2.371 


1.336 


1.239 


1.371 



Table 3: The constants B and Z MS at a — for various levels of improvement. 

8.45 and 8.0, being used in current simulations. The calculations have been performed in 
general covariant gauge, using the symbolic language Mathematica. This gave us complete 
control over the Lorentz and spin structure, the cancellation of infrared divergences, as 
well as the cancellation oil/ a singularities. However, the price is high. In intermediate 
steps we had to deal with O(10 5 ) terms due to the complexity of the gauge field action. 

To improve the convergence of the perturbative series and to get rid of lattice artefacts, 
we have applied tadpole improvement to our results. This was done in two stages. In the 
first stage we improved the gluon propagator, while in the second stage we improved both 
gluon and quark propagators. 

Results at other j3 values, p parameters (also including other gauge field actions with 
up to six links) can be provided on request. 
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